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CN ■ Abstract 

>, 
^O , It is shown that for each finite number N of Dirac measures 6s„ supported 

^ ' at points s„ G M^ with given amplitudes a„ G I^\{0} there exists a unique 

real- valued function u G C'''"^(M^), vanishing at infinity, which distributionally 

r-^ , solves the quasi-linear elliptic partial differential equation of divergence form 

O ; -V • (Vn/Vl-|Vu|2) = 47rJ2n=ianSs„. Moreover, u G C^ {R^\{sn}^=,) ■ 

The result can be interpreted in at least two ways: (a) for any number A^ 
of point charges of arbitrary magnitude and sign at prescribed locations s„ 
in three-dimensional Euclidean space there exists a unique electrostatic field 

^ ' which satisfies the Maxwell-Born-Infeld field equations smoothly away from 

the point charges and vanishes as |s| — )• oo; (b) for any number N of integral 
mean curvatures assigned to locations s„ G M^ C M^'^ there exists a unique 
asymptotically flat, almost everywhere space- like maximal slice with point 
defects of Minkowski spacetime M^'^, having lightcone singularities over the Sn 
but being smooth otherwise, and whose height function vanishes as \s\ — )• oo. 
No struts between the point singularities ever occur. 
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1 Introduction 

In this paper we will prove the existence of unique, essentially smooth distributional 
solutions to the quasi-linear elliptic partial differential problem of divergence form 

W-^=I^l== + AnJ2^rAAs) = for s G M^ (1) 

u{s) — !• as |s| — )■ oo; (2) 

here, 6s^ is the unit Dirac measure supported at s„ G M'^, and the a„ G IR\{0} are 
amplitudes. More precisely, we will prove the following theorem: 

Theorem 1.1. For any finite sets {sn}n=i C M^ and {an}n=i C IR\{0} there exists a 
unique real function u G C°'^(M^) which solves (Qp, ^ in the sense of distributions. 
Furthermore, |Vn(s)| < 1 for s G ]R^\{s„}^^;^, and lim^^^n |Vm(s)| = 1 for each Sn- 
Thus, u G C"^(M3\{s„}Af ^ 



n=l) 



Remark 1.2. Evidently our theorem allows that some of the support points for the 
Dirac measures coincide; however, any such situation is identical to a reformulation 
of the problem with fewer but distinct points, with a re-assignment of amplitude 
values. Thus, without loss of generality we will henceforth assume that all the Sn are 
distinct. The amplitudes may or may not be distinct, though. 

Our result has applications in physics and geometry. It governs objects as diverse 
as, on the one hand, the classical electrostatic fields of the Maxwell-Born-Infeld 
field theory |BoInl934l |Pryl935b[ IGibl998l IKie2004aj . and maximal spacelike hy- 



persurfaces with lightcone defects in the Minkowski spacetime |BaSil98^ IEckl986"| 
IKlMil993| Klyl995 Kly2003| on the other. Applications are discussed in section 4. 



Curiously enough, given the attention that these areas of research have received 
in the literature, the existence of solutions to ([1]), (|2]) as ascertained in Theorem ll.il 
has been an unsettled problem. Of course, there is the explicit solution of ([1]), ([2]) for 
N = 1 found by Born |Borl933j and elaborated on further in |Borl934"t IBoInl93^ 
IBaSil982| IEckl986| IGibl998J : it is well-defined for any value of its amplitude a. 
There is also a semi-explicit solution of ([T]) (which violates ([2]), though) for A^ = oo 
found by Hoppe |Hopl994l and further elaborated on in |Hopl995"t IGibl998j ; it has 



positive and negative amplitude Dirac sources of magnitude \a\ arranged in a cubic 
lattice and exists for arbitrary \a\. However, to the best of our knowledge, generic 
existence theorems for solutions to ([T]) have so far been established only: (a) in 
|KlMil993| Klyl995| under smallness condition^j for the a„ when ([T]) is restricted to 



(bounded or unbounded) domains with boundary with Dirichlet data replacing (|2]) 
(b) in |Kly2003| for arbitrary a„ but with (|2]) replaced by prescribing u{sn) = Ur. 



^We note that the review of Klyachin's paper [Klyl995! in Mathematical Reviews incorrectly 
claims that existence were proved for arbitrary amplitudes. 



restricted by the bounds \un — Um\ < \sn — Sm\ for 1 < n < m < N — in this case 
02]) is generically violated, and it doesn't fohow from the proof in |Kly2003] whether 
(|2]) can hold for some particular choices of {un\n=i '^ ^ ^^^ {0'n}n=i '^ J^\{0}, 
given {sn}n=i C M.^. Our Theorem 11.11 does not follow from adapting the proofs 
in |KlMil993] . |Klyl995| , or |Kly2003| . In fact, our arguments also extend to the 
Dirichlet problem in domains with boundary, as will become clear from our proof. 

As do their proofs of their theorems in |KlMil993j . |Klyl995| , and |Kly2003| , 
our proof of Theorem 11.11 makes convenient use of the results by Bartnik and Simon 
|BaSil982] . Explicitly, in |BaSil982] Bartnik and Simon prove a number of results for 
the Dirichlet problem of ([1]) in bounded domains (with almost arbitrarily irregular 
boundary!), and they also outline how to pass to unbounded domains using barrier 
functions as in |Trel982] . From their results one can extract the following theorem: 

Theorem 1.3. (Essentially Bartnik-Simon.) For any finite set {sn}n=i ^ ^^ ^/ 
distinct points and any finite set {un}n=i C M, restricted by the bounds 

\Un — Um\ < \Sn — Sm\ for 1 < U < TTl < N, (3) 

there exists a unique real function u G C°'^(M'^) which weakly solves 

V ■ , ^""l"] ,„ =0 for se R\{s4Li, (4) 

^1- \Vu{s)\^ 

u{s) —> Un as S — )■ Sn, (5) 

u{s) — > as \s\ — )■ oo. (6) 

Furthermore, \Vu{s)\ < 1 for s e R^\{sn}n=i, and u G C^(M2\{s„}^^i). 



Theorem 11.31 basically reduces the proof of our Theorem 11.11 to variational argu- 
ments which show that for each set of amplitudes {an}n=i C IR\{0} associated with 
the points {sn}n=i C M.^ there exists a unique distributional C°'^(]R'^) solution u{s) 
of ([I]), ([2]) for which |u(s„) — u{sm)\ < |s„ — s^l for 1 < n < m < A^. The claim that 
lim|5_s„|^o |Vm(s)| = 1 then follows from: 

Theorem 1.4. (Rephrasing of Theorem 1.5 in \Eckl986^ ). Letu{s) be as in Theorem 
\1.3[ Then either u{s) can be analytically continued into s„ or 

lim |Vu(s)| = 1. 

Thus, genuine singularities of u{s) are lightcone singularities. 

In section 2 we formulate our variational approach to ([1]) , ([2]) and prove existence 
of a unique optimizer u E Cq' (M^) with IVmI < 1. In section 3 we show with a 
dual variational argument that |u(s„) — u{sm)\ < \sn — Sm\ ioi 1 < n < m < N. 
Afterwards, in section 4, we discuss applications to physics and geometry. In section 
5 we list a few straightforward extensions of our main theorem, only indicating their 
proofs. In section 6 we conclude with a list of desirable extensions. 



2 The variational approach 

In this section we prove: 

Proposition 2.1. There exists a unique u G Cq (E?) fl {v : |Vf | < 1} for which 



l<n<N 



holds for all ip G C^(M^), where Vm denotes weak derivative and where d'^s zs three- 
dimensional Lebesgue measure. 

Thus, u e Co'^(M3) n {?; : \Vv\ < 1} is a distributional solution of (^, (jS]). 

2.1 Preliminary considerations 

Equation ([1]) is the formal Euler-Lagrange equation for the variational principle 

3'iv)= I (l-v/l-|Vt;(s)|2-47rt;(s) J^ a„(5,„(s))d=^s ^ min (8) 

over a suitable set of functions f , and ([Tj) is its weak version. In particular, if C^(]R^)n 
C^(M^\{s„}^^]^) denotes the Banach space of bounded continuous real functions on 
M'^ which have a bounded continuous derivative on the indicated punctured domain, 
equipped with their usual norm, then J is well-defined for those v G C°(]R^) fl 
C^(M^\{s„}^^;^) which satisfy |Vf | < 1 on ]R^\{s„}^^]^ and which vanish sufficiently 
fast as |s| -^ oo; in particular, |Vf(s)| = 0(|s|~^) is fast enough. Eventually, in 
section 3, we will show that J does take its minimizer on this set of functions. 
However, since the indicated Banach spaces are not convenient spaces to work with, 
here we shall characterize ([8]) as upper limit of a sequence of variational functionals 
which are defined on larger, more convenient spaces of functions. The minimizer 
of 5" will be obtained as the limit of a family of minimizers of these approximating 
variational principles. In particular, we show that the minimizer solves ([7]). 

2.2 A monotone family of variational principles 

For X > we define the extended real- valued function 



/ 1 - \/l-x for X G [0, 1] , 

(^ oo for a; > 1 

The i^-th Taylor polynomial of / about a; = 0, given by 

K 

Tay,,[/](x|0) = 5^/('=)(0)x^ (10) 

A:=0 



has Taylor coefficients 

for A; = 

/('=)(0)=<' \ for /.= ! , (11) 



(2fc-3)! 



for k > 1 



{2k)\\ 

SO that Tay[/](a;|0) = {Tay;^[/](x|0)}^^^, the Maclaurin series of /(x), viz. 

Tay[/](x|0) = ia: + f;^^|^x^ (12) 

is a pointwise strictly increasing sequence of strictly convex, strictly increasing func- 
tions of a; > which vanish at a; = 0. The series flT2|) converges absolutely to f{x) for 
a; G [0, 1] but diverges for x > 1; since all coefficients are positive, Tay^[/](x|0) /^ oo 
for a; > 1, so we are entitled to say that Tay[/](a;|0) actually converges to the ex- 
tended real- valued function f{x) for all a: > 0. In the following, for brevity we shall 
simply write fx^x) for Tay^[/](x|0). 

With the help of the Taylor polynomials we now define the family of functionals 

'^k{v)= [ (fK{\Vv\^)-^nv{s) V aj,„{s))dh, (13) 

which for K > 2 are well-defineco on ni<fc<i<:^o' (I^^)- For the |Vf| term this 
is seen right away from the definition of W^' (M^) as the closure of the compactly 
supported C°° functions on M'^ w.r.t. ||f ||p^i,2fe,jg3> = /j^g |Vv(s)p^d^s, so that for 



'0 

V G ^x<k<K^l^'^^^) we have /i^dVt;^) G ^^(M^). To see that also the source 

term in flT3|) is well-defined we note that Wg' (^^^) '^ ^lo'c (^^) ^° t'^^^ ^^ ""^^ apply 
Sobolev's original embedding theorem, according to which for any ball 5 C M^ we 
hav(i W^^'^^{B) ^ C^{B) whenever fc > 2, and conclude that 5^^ G W^^'^'^^^\^^) 
for all k > 2. This estabhshes that 3^k is well-defined on ni<fc<A'^o' (^^) ^^^ 
K > 2, and any finite N. Incidentally, since elements of Wq' (M^) tend to zero at 
spatial infinitjQ a.e., we see that the v G (~]i<k<K^o' (^^) with K > 2 are in fact 
equivalent to a subset of the bounded continuous functions on M^ which vanish at 
spatial infinity. 

Let =e/ = ni<fc<oo^o' (I^^)n{f : |Vf| < 1} be the set of admissible (equivalence 

classes of) functions for 3". Since s^ C ni<A,<i^^o' (^^) ^^^ all K e N, and 
since fxix) /^ f{x) for all x > 0, monotone convergence now yields that 3^{v) = 
lim/i'^oo 3^Kiv) for all f G =2/. 



^The functional J'l is not well-defined on the "canonical" domain of the Dirichlet integral, which 
is Wq ' (M^), for which reason we don't have any use for J'k when K = 1. 

^By the Sobolev-Morrey embedding theorem we even have W^'^''{B) '-> C'^'^~^^^''{B) for fc > 2 
and any B CR^. 

''This is not true for ah elements of Wq''^''{R^) when k>2 



2.3 Existence of a family of unique minimizers 

For each K > 2, the functional 3^k is clearly convex over ni<fc<_ft:^o' (^^)- More- 
over, 3^K is bounded below and coercive w.r.t. the topology of ni<fc<E'^o' (^^) 
whenever K > 2. To see this we have to estimate the source term in 3^k- 
We rewrite 3^k as 

K 

^k{v) = '^Ck\\vf^i,2k^^3^-4:n ^ anv{sn), (14) 

fc=l l<n<Af 

where Ck = f^''\0) > 0. Now {sn}n=i is given, so there exists an open ball i? C M^ 
such that {sn}n=i C B. Since the restriction of any v G Wq' (M^) to B is in 
W^'^{B), and since W^''^{B) C W^'^{B) (though no embedding, clearly), we can 
apply the Sobolev embedding theorem in the form W^'^{B) "^^ C^{B), and using 
|a„| < maxi<„<Ar |a„| < oo, for all i^ > 2 we obtain the estimate 

K 

^k{v) > Y^ Ck\\v f^i,2k ^^3) - 47riVA||t;||^^i,4(s) (15) 

fc=i 

where A is a positive constant. Now ||'y||^;(^i,4/^-) = ||'y||j^4(-m + ||Vf ||j;^4(-m, so that 

1 1 '^ 1 1 ly 1.4(B) — il'^llL4(m + l|Vf 11^4(3)- Next, since v G Wq' (M^) for K >2, we have the 
nontrivial estimate ||Vf||^4(£) < ihllvi/i'-ims) < oo. Furthermore, by Holder's inequal- 
ity, II^'IIl4(b) < \B\^/^^\\v\\^,(^j^-^, and the special case Wq'^(M.^) ^ L^{R^) of Sobolev's 
embedding theorem then yields the estimate ||'y||j^4(m < |-B|^/^^S'||f ||p^i,2.jjj3x < oo, 
where 5* > is the sharp Sobolev constant. And so, for K >2 we find 

K 

^^Kiv) > J2 Cfcll^llS^J-2fc{R3) - 47rA^ (^'II^IU,}-2(R3) + Mv\\w^:^^K3)) , (16) 

fc=i 

where A' is another positive constant. This lower estimate for 3^^ is manifestly 
bounded below and coercive on ni<fc<i^^o' (^^) whenever K > 2. Therefore, 
for each K > 2 the functional 3^k takes on a unique minimum for some vk G 
f]i<k<K^o'^\^')- We set 3^KivK) = Fj,. 

2.4 Weak convergence of the family of minimizers 

Since F^ = 3^k{vk) > 3^k'{vk) > 3^k'{vk') = F^' whenever K > K', the minimum 
values Fk of the family of variational functionals form a strictly monotonic increasing 
sequence. And since fxix) < f{x) for all K when a; > 0, this sequence {Fk}k=2 is 
bounded above by 3^{v), where v G C"^'^(]R^) is the following convenient trial function: 
let 2r := min{|sfc — si\}i<k<i<N, then v is defined by 

^ f •. _ j s^g^ (an) {r - \sn - s\) for s e Br{Sn) ..,„s 

''^'^"l for 5GM3\Ui<n<^i?r(s„) ^''^ 

6 



One readily calculates that 

3^{v) = N (\Br\ - Anr'mY (18) 

where 

l<n<7V 

Thus, linij^^oo Fk =: F < inf^ J(t;) < 3^({)) exists, and Fk < F for all K>2. 

As a corollary, since S^k'^vk) < 3^k{vk) when K > K' , we have that 5'i<:'(t'ii:) < -F 
whenever K > K' . By coercivity, for any fixed K' > 2 there now exists a positive 
constant C such that ||fii-||^i,2K'.jg3. < CF for all i^ > K'. Now, since Wg' (M^) 

is a separable, reflexive Banach space for all 1 < i^' < oo, the closed ball |w : 
||'y||^i,2K'.jg3s < CF} is weakly compact. Therefore the sequence {vk}k=2 contains 

a weakly convergent subsequence in (VFg' ^ ^o' )(^^) ^^ each K' > 1. By a 
diagonal argument we can pick the subsequence so that its weak limit in each {Wq ' fl 
Wo''^^'){R^) is one and the same Woo G r\i<k<K' W^'^^R^) for all 1 < K' < oc, hence 

^ooeni<.<ooW^o'''(K')- 

Now, functions in ni<fc<oo ^o ' (-l^^) ^^^ ^^^ necessarily in PVq '°°, but the uniform 
(in K) upper bound F on the 3^k{vk) guarantees that the weak limit Woo of the vk 
is actually in Wq'°°; indeed, we even have |Vfoo| < 1 a-e. For assume to the contrary 
that |Vfoo| ^ 1 a.e. Then there exists an fi C M^ with |r2| > such that |Vfoo| ^ 
l + 2e a.e. in Q. But then there exists a K such that | Vf/^l > 1 + e a.e. in Q whenever 
K> K. And then we have S^Kivx) > \^\ T.k=K+i ^^^(1 + ^)^'' / oo as iT ^ oo, 
which is a contradiction to 3^k{vk) < F for all K. Thus, not only is Voo G Wq'°°, also 
iVwool < 1 a.e. in M^ as claimed. So i^oo G ni<fc<oo Wq'^''{R^) n {t; : |Vw| <!} = £/. 

2.5 The limit of the minimizers of the 3^k minimizes 3^ 

We have just proved that foo G i^/. We now show that 3'(t'oo) = min^jg^ 5'('^)- 

Since j/ C {W^'^ D W^'^^)(R^) for each K > 1, g^A'(^) is well-defined for each 

V & £/ and each K G N; in particular, 3^k{vqo) is well-defined for each K E N. 

Moreover, limK^oo3^K{v) = 3^{v) for each v G ^, by monotone convergence. 

Now, by the monotonicity of the sequence of Taylor polynomials {/ii'dVf |)}^^-^, 

we have 5'(foo) > ^^^(^oo) > 3^k{vk) = Fk for all K > 1. By taking the limit as 

i^ — 7- oo, we obtain 3'(foo) > lim_ft:^oo3'E'(^A') = -F- 

On the other hand, recalling flT^ . we see that each 3^k is obviously weakly 

lower semi-continuous, so we have 3^k'{voc) < ^iT^K^oo3^K'{vK)- Recalling now that 

^^k'^vk) < F whenever K > K', we obtain 3^k'{voo) < F for all K' > 1. Taking the 

limit K' -^ oo and recalling that lim^^/^oo ^k'{v) = J^{v) for each v G ^, we obtain 

that 3'(woo) < -F- In total, we have shown that 5'(t'oo) = F. 

7 



It remains to show that there is no {; G =2/ for which 3^{v) < F. But this is really 
easy. For assume to the contrary that there were such a v with 3^{v) = F — e. Then 
J^xiv) < F — e for all i^ > 1, which contradicts the fact that for any e we can find a 
K{e) such that min^£(vi/i,2niyi>2K){K3) 3^k{v) = 3^k{vk) > F — e whenever K > K{e). 

This proves that 5" takes its minimum at foo G sz/, and the minimum equals F. 
Moreover, by convexity, the minimizer is unique. 

2.6 The minimizer of 3^ weakly satisfies the Euler-Lagrange equation 

We cannot yet conclude that the minimizer foo of 3^{v) weakly satisfies the for- 
mal Euler-Lagrange equation ([T]) because for this conclusion we need to know that 
|Vfoo| < 1 a-e. and so far we only know that |Vfoo| < 1 a.e. We now show that 
|Vfoo| < 1 a.e., which implies that foo weakly satisfies ([1]), i.e. ([7]). 

Let i^crit = ne>o {^ '■ l^'^ool > 1 ^ c}- Notc that ficvit contains all points s* at 
which I Vf 00(5*) I = 1 as well as all points s^ for which ess-lims_j.s^ |Vfoo(s)| = 1 
without necessarily having Vvoo{s) itself defined at s = 5=,,. Clearly, ficnt has finite 
Lebesgue measure, |f2crit| < C)0, for f oo G ^ implies that iVfool-s)] — ?■ as |s| — )■ oo. 
We now show that |f2crit| = 0. 

For this purpose we assume to the contrary that |f2crit| > 0. Then the variation 
of 3'(f ) about foo gives to the leading order (i.e. power 1/2) in ip, 

^^'/'^boo](V'):=- / y/-2Vvoo{s)-Vil:{s)dh, (20) 

where ip ^ C^(R'^) is any test function satisfying Vfoo(s) ■ Vipi^s) < a.e. on ficrit- 
Note that 3^*-^''^'' [f oo] (V") is homogeneous of fractional degree 1/2 in -ip; hence, this 
nonlinear term — if nonzero — will in general dominate the usual linear terms in 
tjj, indeed. Moreover, whenever 3^^^^'^^[voo]ifp) 7^ 0, we manifestly have 



y/-Vv^{s)-V^{s)d's<0. (21) 

But for foo to minimize 5" over jz/ we must have 3^^^^'^^[voo]{i') > for all 'ip E C^(M.^) 
satisfying Vfoo(s) ■ Vip{s) < on i^crit a.e. This together with ( 121]) implies that 
3^(1/2) [^^](^) = for all ^ e C^{R^) satisfying Vv^{s) ■ Vip{s) < a.e. on ficrit- 
But this is only possible if |f2crit| = 0, as claimed. 

Remark 2.2. Our argument above does not show that i^crit = {sn}n=i- 

The result |i7crit| = means that |Vuoo| < 1 a.e., and this already implies that 
the variation of 5'(f ) about Voo to leading order (i.e. power 1) in ip now reads 

3'^'\v^m= [ (ViPis) ■ ^=I^d^l== - AmPis) Y^ ^nSsAs))d's. (22) 



Since J'^^^[voo]{ip) is linear in ip, Vao can minimize 5" over £/ only if S'^^^[voo]{'ip) = 
for all Ip, which is precisely ([7]). Thus the Euler-Lagrange equation ([T]) is satisfied by 
foo in the weak sense, as claimed. The proof of Proposition [2]T] is complete. Q.E.D. 

Remark 2.3. We close this section with the remark that alternate, nonvariational 
routes to Proposition \2.1\ are conceivable. In particular, the Dirac sources can he 
mollified with compactly supported C°° functions, and the asymptotic vanishing of 
u{s) as \s\ — )■ oo replaced by 0-Dirichlet conditions on OBr, where R is a large ball 
containing the supports of all mollifiers of the Dirac sources. For this situation the 
theorems in \BaSil98^ guarantee a classical solution to the so mollified (QP- As 
pointed out by one of the referees. Lemma 2.1 in \BaSil982^ and elliptic regularity 
theory should now yield uniform Lipschitz bounds on the solution u away from the 
eventual locations of the Dirac sources when the mollifiers are removed, and the 
proof of their Lemma 3.1 shows that the limit function solves ^ restricted to Bfi. 
Subsequently one can let R ^ oo by invoking Treiberg's barrier function arguments. 

3 Bootstrapping regularity 

In this section we bootstrap the regularity of the minimizer i;oo = u oi 3^{v) to the 
level which guarantees satisfaction of Theorem II. 1[ 

3.1 Bootstrapping the regularity of m away from ficrit 

By our Proposition 12.11 the unique distributional solution to ([T]), ([2]) obtained by 
minimizing 3^ in ^ takes values «„ = u{s„) at the s„ which satisfy the inequalities 
l^n — Um\ < |s„ — Sm| for all 1 < n < m < A^. Hence we can invoke Corollary 4.2 to 
Theorem 4.1 of [BaSil982] to extract the following proposition for our setting. 

Proposition 3.1. Let u{s) G C°'^(M'^) be the unique distributional solution to (|7]j, 
^ which minimizes 5" in s^ . Then u G C"^(R^\r2crit)- Moreover, Vtcrw (the singular 
set K for u in WaSil98^ ) is a subgraph of Xn = 'X.(^{sn}n=i) C M^, the complete 
graph whose vertices are the set {sn}n=i- Furthermore, let En,m C %n denote the 
edge of %n with endpoints s„ and Sm- Then En^m C r^cnt if and only if \un — Um\ = 
\sn — Sm\, and in that case we have u{tSn + (1 — t)sm) = tUn + (1 — t)um for t G [0, 1]. 



3.2 Proof that Vtcrit = {■§„} 



N 
n=l 



We recall that any distributional solution G W^'"^ of ([T]) satisfies the weak maximum 
principle, Theorem 8.1 in |GiTrl98^ . Therefore foo(s) = u{s) has a local maximum 
at Sn whenever a^ > and a local minimum when a„ < 0, and no extremum in 
]R^\{s„}^^]^. This together with Proposition 13. II right away gives us: 

Corollary 3.2. Let ana^ > 0. Then Enm\{sn, Sm} 't- ^crit- 



Thus, the only potentially critical edges -E'n^m are those whose end points s„ and 
Sm sport amplitudes a„ and a^ of different sign. To show that also those edges, 
save their endpoints, are not critical requires a different argument. We shall invoke 
a convex duality argument which rules out all the edges, save their endpoints, from 
the critical set. 

Proposition 3.3. For all 1 < n < m < N we have that Enm\{sn, Sm} ^ ^crit- 

Proof: Since C^(M^) is dense in Cq' (M'^), we can substitute Voo = u for ■?/' in ([7]) 
and, for the solution m of ([7]), obtain the identity 

0= / ( ;'^"!!^'\„ -4^^(^) E ^nSsAs))d's. (23) 

A simple rewriting of f l23|) yields 

3^(u) = [ (l- , ^ )d^g. (24) 

Defining 

Uis) ^ ,.-^;W ^, . (25) 



where m is still the solution of ([7]), another elementary rewriting of f l24p yields that 
J{u) = -g(f/), where 

g(l^) = [ [VI + I^P- lid's (26) 



is well-defined for any vector field V for which |l^| G Lj'qj,(M^) fl L^(M^\i?ij), where 
-B/j is some ball of large radius R. Note that 9(y) is related to 3^(f ) by a Legendre- 
Fenchel transform, viz. 

9(V) = max /" f[v/l-|Vt;P- ll - ^ ■ Vv) dh; (27) 

the dual variables of the transformation are Vf <-> V. Thus, 9{V) is strictly convex 
in V. But we have seen that also 3^{v) is strictly convex for v G ^, so 3^(f ) — or 
rather the source-free part of 3^{v) — is given as a Legendre-Fenchel transform of 
S{V). As a result, we can also obtain the minimum of 3^{v) and its minimizer Voo = u 
in terms of a constrained minimum principle for 9{U). Explicitly, 

r ^ ^ 

9{U) = min 9iV) V-V = AttY: aJs^; \V\ G Aoc(K') n L2(R3\5«) ; (28) 

•- n=l ■' 

in (128!) it is understood that V ■ V^ is well-defined in the sense of distributions and 
that R is big enough so that {sn}n=i C Br. 
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Next, we define the almost everywhere harmonic field 

N I 

V^/.(s) = -E«nV^^. (29) 

n=l I"' "'ral 

Note that 14(s) G Aoc(^^) ^ L^{'^^\Br) whenever {sn}^=i C Br. Furthermore, 

AT 

V-Vh = A7T^aJsr.. (30) 

ra=l 

So V/i is in the set of admissible vector fields for our variational principle ( 128|) . 

We are now ready for our main argument. Namely, suppose that for some n,m 
the edge En^m C ficrit- Without loss of generality we may assume that no other 
Sk lies on En,m- Then lim^^£;„^ |Vm(s)| = 1, and so hm^.^^;^^ \U{s)\ = oo. But 
since u is analytic away from f^crit, so is U, hence we conclude that there is some 
tubular neighborhood of En^m in which |t/(s)| > |V/i(s)|. Since s„ ^ Sm we can 
intersect our tubular neighborhood of En^m with two small balls centered at s„ and 
Smi respectively, and delete the intersection domain from it. Denote the resulting 
truncated tubular neighborhood by -E"^; it is a bounded open set. Mollify its 
boundary (9i?° ^ a little bit to obtain an open corridor C° ^; it needs to have a finite 
distance from any Sk- Now let K(s) be given by U{s) ioi s ^ E^^U C^mi ^^d by 
K(s) = Vh{s) for s G -E'°^\C°„. We need to connect these fields smoothly across 
C°^, but this is easy. Since away from {sn\n=i the fields U and Vh are divergence- 
free, we can represent each field as the curl of some vector field. We can choose a 
C°° deformation of one such vector field into the other across the transition region 
C° ^, and in C° ^ we define K to be the curl of this deformed field. Thus V^ is in 
the set of admissible vector fields for our variational principle. But then we have 
S(K) < S(f^), which is a contradition to our variational principle (125]1 . 

Thus \U{s)\ < oo for s G M\{s„}^^=i, and therefore ficrit = {sn}n=i- Q-E.D. 

Remark 3.4. We close this section with the remark that our convex duality argument 
can also he adapted to show that |Vvoo(s)| < 1 away from {sn}n=i without invoking 
Proposition \3 . 1\ But then a Nash-Moser argument has to be supplied to bootstrap the 
regularity of Vqo from Lipschitz continuity to real analyticity in M.^\{sn}n=i- 

4 Applications to Geometry and Physics 

4.1 Spacetime interpretation of Theorem [TTTl 

A smooth space-like hypersurface S in Minkowski spacetime M^ = M x M'^ is a three- 
dimensional simply connected subset of M^ with a time-like normal vector at every 
point in E. Thus S = {zu G M^ : 7{to) = 0} is the boundary of the zero level set 
of a differentiable function T : M^ — > M with ran(T) = M and dT(ct7) time-like, i.e. 
g~^(dT(tz7), d'J(ci7)) < for all zu G M^; here d is E. Cartan's exterior derivative on 
M'^ and g the Minkowski metric with signature +2, a 2-covariant tensor acting on 
T(M^) X T(M^), where T(M^) is the tangent bundle of M^. Topological^, S ~ R^. 
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Since any such hypersurface is a graph over M^, without loss of generahty we may 
assume that the generating function T is of the form T{w) = t — c^^S{s), where 
zu = {ct, s) defines a Lorentz frame, where t is time and s is a vector in Euchdean 
space M3. ThenE = {(ct,s) : t = c-^S{s)}. Since g-\dT (tu) , dT (w)) = -l + |V5|^ 
and since S is space-hke, we need to have 1 — \'VS\ > everywhere. 

For those S which are asymptotically flat, more precisely if S x Eg with Sq = M^, 
the volume difference Avol(E|Eo) of E versus Eq is well-defined. After at most a 
Lorentz transformation we can choose Eq = {{ct, s) : t = 0} = M^, in which case 



Avol(E|Eo) = / [yl- I V5r - 1 J (Ts. (31) 

Note that Avol(E|Eo) < 0. A hypersurface E is called maximal if any compact 
variation leads to a decrease of volume. In particular, Eq = M^ is a maximal entire 
spacelike hypersurface in M'*. By a Bernstein theorem of Cheng and Yau |ChYal976] . 
any entire space-like maximal hypersurface in M^ is fiat; see also |Yan2000] . Thus, 
to be interesting a maximal hypersurface cannot be entirely space-like but at best 
only space-like almost everywhere. If in particular E has isolated defects then by 
Ecker's theorem these are lightcone singularities, i.e. isolated points in E where the 
normal vector touches the lightcone. Any such almost-everywhere space-like maximal 
hypersurface with point defects in M'' is the graph E = {{ct,s) : t = c~^S{s)} of a 
function S{s) satisfying 1 — \VS\ > away from the defects, such that 1 — |VS'| 
extends continuously into the defects, where it vanishes. 

Prescribing the locations Sk G M^ of the lightcone singularities does not uniquely 
determine an asymptotically flat maximal hypersurface with defects. In addition, 
the particular asymptotically linear behavior of S{s), and also the integral mean 
curvatures /ife G IR\{0} which are associated with each lighcone singularity of the 
hypersurface have to be prescribed. Maximizing Avol(E|Eo) for such a hypersurface 
with lightcone singularities of prescribed integral mean curvatures is a variational 
problem with constraints. The Euler-Lagrange equation for this problem readqj 

-V ■ , = 3 V fikSs, . (32) 

Identifying S = u and Hk = (47r/3)afc yields ([1]). 

In this spacetime interpretation our Theorem 11.11 becomes: 

Corollary 4.1. For any set {sfcj^i C M'^ of distinct points and any set of integral 
mean curvatures {fik}k=i C IR\{0} assigned to these points, there exists a unique 
asymptotically fiat hypersurface E = {(ct, s) : t = c~^S{s)} with S E Cq' (R^,R) fl 
C"^(R2\{sfc}f^i,]R) solving ^; moreover, \VS{s)\ -^ 1 as s ^ Sk- Thus E is an 
almost everywhere space-like maximal hypersurface, having N lightcone singularities 
with prescribed integral mean curvatures fik located at the Sk- 



^We follow the convention of |GiTrl983] which differs from that in [BaSil982| by the factor 3. 
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4.2 Electrostatic interpretation of Theorem 11.11 

In classical electromagnetic field theory, the Coulomb law states that an electric 
(point-) charge "density" in M.^ is the source of the electric displacement field D, 

N 

V ■ D = Any^Zkds^ Coulomb's law , (33) 

fc=i 

witlij Zk E Z\{0}, while Faraday's law says that an electrostatic field E is curl-free, 

V X E = Faraday's law (stationary) . (34) 

These two laws need to be complemented by a law relate E and D, for which Max 
Born |Borl933J proposed 

D = — — Born s law (35) 

^l-f3^\E\^ 

with P G (0, oo) (we use the dimensionless notation of |Kie2004a] ). In the limit 
/3 — )■ Born's law fl55]) goes over into Maxwell's law of the "pure aether", D = E. We 
remark that (13^ is the electrostatic limit of both, the electromagnetic law proposed 
by Born |Borl933l IBorl969] and the more elaborate law proposed by Born and Infeld 
|BoInl93^ IBoInl93^ . The latter has received much attention in recent years, see 
the references in |Gibl998l IKie2004al IKie2004bl[Kie2012j . 

Clearly, (|34ll implies that E = —VA for some scalar potential field A. Inserting 
this representation for E into ( l35l) . which in turn is inserted in (l33ll . we obtain 



V4 Jl 

-V ■ , = inyzkS,, . (36) 

Multiplying (15^ by /3^ and identifying [3"^ A = u and P'^z^ = a^ we retrieve ([T]). 
In this electrostatic interpretation our Theorem 11.11 yields: 



Corollary 4.2. For any finite number N of point charges {z^}^^-^ C Z\{0} located at 
distinct points {sk}^^i C M^, there exists a unique electrostatic field E inM^^\{sk}^^i 
which solves (E^j; ^i34\ ), (E^) cind has finite field energy^ 

1 a 



«»-(D) = -^y^_(VlH-flW-ljd^«. (37) 

The solution E G C"^(]R'^\{sfc}^]^,]R^), hut it cannot he continuously extended into 
the Sfe. It is hounded, with (3'^\E{s)\ — )■ 1 for s — )■ s„, and it vanishes for \s\ -^ oo. 



^Empirically, all nuclear and electron charges are integer multiples of the elementary charge, 
which is unity in our units. 

^Here, a is Sommerfeld's fine structure constant, inherited from the units in |Kie2004a| . 
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Remark 4.3. Presumably inspired by Theorem 4-i and Corollary 4-2 in lBaSil982^ . 
at the beginning of section 4 in \Gibl998 ln Gibbons contemplates the following: "It 



is well known that one can construct explicit multi-black hole solutions held apart 
by struts, the struts being the sites of conical singularities representing distributional 
stresses. One should be able to construct analogous multi-BIon solutions. " (What 
Gibbons calls "multi-BIon" solutions are but electrostatic solutions to the Maxwell- 
Born-Infeld equations with many point charge sources. In particular, Bom's solution, 
for Gibbons, is "the Blon.") Our Theorem \l.l\ and its Corollary \4 . S\ show that struts 
between the point charges do not occur in the electrostatic solutions to the Maxwell- 
Born-Infeld field equations with point charge sources. 

5 Extensions of our results 

The geometric interpretation of u as time function of a maximal almost-everywhere 
space-like hypersurface with lightcone singularities in Minkowski spacetime allows 
one to exploit the Poincare group of M^ to generate solutions m to ([1]) with different 
linear asymptotics at space-like infinity than ([2]). Since there is a unique Poincare 
transformation for the transition from asymptotically vanishing to non-zero asymp- 
totically linear conditions, we can therefore conclude: 

Corollary 5.1. For any finite sets {sn}n=i C M'^ and {a„}^^^ C IR\{0} and a vector 
e ^M.'^ of magnitude \e\ < 1 there exists a unique real function u G C'^'"'^(]R'^) satisfying 

u{s) — e ■ s — !■ as |s| — )■ oo, (38) 

and solving (QP in the sense of distributions. Furthermore, |V'u(s)| < 1 for s G 
^^\{sn}^=i, and lim,^,„ |Vm(s)| = 1 for each Sn. Thus, u G C'^(M3\{s„}^=i). 

The equivalence between the mathematical theories of maximal space-like hyper- 
surfaces with point defects in Minkowski spacetime on the one hand and the electro- 
static Maxwell-Born-Infeld potentials generated by point charge sources on the other 
allows us furthermore to re-interpret these asymptotically nontrivially linear hyper- 
surfaces as electrostatic solutions with asymptotically (at spacelike infinity) constant 
electric fields. It is worth stressing that, in the notation of our previous subsection, 
one thus interprets (/3^A, s), rather than the spacetime point {ct,s), as Minkowski 
four-vector to generate new solutions by Poincare transformations. This "hidden 
Poincare symmetry" seems to have been exploited first by Gibbons, see sections 3.3 



^In the same paragraph in jGibl998] the results of Bartnik and Simon |BaSil982) (Gibbons' 
reference |Barl987] ) for the weak solvabiUty of the Dirichlet problem are somewhat misquoted: the 
necessary condition (in our notation) |u„ — Um\ ^ |sn — Sm| on the Dirichlet data at the distinct 
points Sn and Sm, which in Gibbons' notation would have to read |<i>° — $''| < \xa — Xb|, is missing. 
Whenever |$° — $''1 — [x^ — X{,|, then a strut between Xa and Xf, does occur. 
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and 3.7 of |Gibl998] : in particular, in section 3.7 Gibbons transforms Born's solu- 
tion into an electrostatic solution with a single point charge and an asymptotically 
constant electric field whose magnitude is below Born's critical field strength. 

Lastly, as already announced in the introduction, there is an analogue of our 
Theorem 11.11 for the Dirichlet problem in bounded domains with nice boundary. 
This is not directly a corollary of our proof, yet its proof follows by a straightforward 
adaptation of our proof to the Dirichlet problem. Thus we claim: 

Theorem 5.2. Let Vt dM? he a hounded Lipschitz domain containing the finite point 
set {sn}n=i C M^. Then for any 4> '■ dQ — t- M satisfying \(p{s) — (p{s')\ < \s — s'\ for 
s 7^ s' , and any set {an}n=i C ]R\{0}, there exists a unique real u G C^'^iVt) solving 

V ■ , 2^\ ^,^ + 4vr 5^ a„(5.„ (s) = for s e Q, (39) 

u{s) = (p{s) for seOn (40) 

in the sense of distributions. Furthermore, |Vm(s)| < 1 for s G fi\{s„}^^i, and 
lim^.^^,^ |Vm(s)| = 1 for each s„. Thus, u G C"^(r2\{s„} 



n=l) 



Remark 5.3. Bernd Kawohl kindly explained to me that for such a hounded domain 
the detour via the 3^k{i^) should he unnecessary to minimize the convex functional 
3^{v) over the convex set {v G Wq'°°{Q) : |Vf| < la.e. in Q}. 

6 Desiderata 

For matters of a quantitative nature it is important to have efficient algorithms to 
actually compute the solutions which in this paper we have proved do exist. For 
Holder-continuous regularizations of the point charge sources such an algorithm has 
been developed in [CaKi2010t IKie2011] , but so far none seems available which would 



cover point charge and other singular sources in R^. The situation is better for 
the lower- dimensional problem in M^, see [Pry 1935a IKobl988| IFer2010J . and it is 



desirable also for the solutions in M.^ to have explicit formulas in terms of, say, 
quadratures and such. For the time being, the variational arguments allow one to 
work with numerical discretizations and to run minimization routines. 

Another question is whether our Minkowski space results extend to certain curved 
Lorentz manifolds, in particular to asymptotically flat Lorentz manifolds |Barl984] . 
If the Lorentz manifold is given (a so-called background spacetime), then the essence 
of the results of [BaSil982] remains true under appropriate conditions, as shown 
by Bartnik in [Barl988] with quite different arguments than those in [BaSil982j . 
Moreover, in [Barl989] Bartnik also extended Ecker's singularity theorem to certain 
Lorentz manifolds. For those Lorentz manifolds for which the analogue of the fiat 
spacetime theorems of Bartnik-Simon hold we expect that analogues of our theorems 
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will hold as well. We remark that Bartnik's theorems in |Barl988"| IBarl989] do not 
require the manifold to be asymptotically flat. 

Another question, of prime importance as explained in |Kie2012] . is whether 
the extension of our electrostatic results to a general-relativistic setting is possible 
in which an asymptotically flat Lorentz manifold is to be found by solving Ein- 
stein's fleld equations, with an electrostatic energy(-density)-momentum(-density)- 
stress tensor as curvature source for the metric, along with solving the Maxwell- 
Born-Infeld equations for the electrostatic fleld in the curved spacetime. The prob- 
lem with a single point charge source was treated already by Hoffmann |Hofl933^ 
but only recently with complete rigor, by Tahvildar-Zadeh |TaZa201l] : there the 
reader is also directed to the large literature on the subject. The general sentiment, 
as expressed in the quote from Gibbons at the end of section 4, seems to be that in 
the multi-point-charge problem struts will occur between the point charges; see also 
|Weil996] . We hope to offer a deflnitive answer in the foreseeable future. 
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